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Abstract 

We investigate the electroweak phase transition in the presence of a large CP 
violation in the squark sector of the MSSM. When the CP violation is large, scalar- 
pseudoscalar mixing of the Higgs bosons occurs and a large CP violation in the 
Higgs sector is induced. It, however, weakens first-order phase transition before the 
mixing reaches the maximal. Even when the CP violation in the squark sector is 
not so large that the phase transition is strongly first order, the phase difference 
between the broken and symmetric phase regions grows to 0(1), which leads to 
successful baryogenesis, when the charged Higgs bosons is light. 
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1 Introduction 


The baryon asymmetry of the universe (BAU) is one of the most obvious facts, which has 
been a longstanding problem in astrophysics[l]. To explain the light-element abundances 
within the framework of the standard big-bang nucleosynthesis, it is required that [2] 

^ = (0.21-0.90) X 10“^°. (1) 

It is well known that in order to obtain this asymmetry starting from a baryon-symmetric 
state, three requirements must be satisfied: baryon number violation, C and CP violation, 
and departure from equilibrium[3]. In general, electroweak theories satisfy baryon number 
violation through chiral anomaly and have a possibility to generate the BAU[4]. In the 
minimal standard model(MSM), the main source of CP violation comes from the phase 
(5 km in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. Although this 
phase is able to account for the experimentally observed CP violation in the neutral K- 
mesons and, as recently measured, in the system, it has been shown that it is not 
possible to generate sufficient BAU through (5 km- Furthermore, the strength of the phase 
transition is so weak in the MSM with the Higgs scalar heavier than 115GeV[5, 6] that 
the nniverse is approximately in eqnilibrinm, when baryon nnmber changing process is in 
effective. 

In the context of snpersymmetric (SUSY) extension of the MSM, it has been pointed 
ont that in the presence of a light stop the electroweak phase transition (EWPT) can be 
strong enongh for baryogenesis to take place[7]. Moreover, SUSY models contain many 
complex parameters as new UP-violating sonrces in addition to (5 km; the Higgs bilinear 
term, fi and soft SUSY breaking terms (gangino masses and scalar trilinear conplings) [8]. 
Besides these complex parameters, the relative phase of the expectation valnes of the two 
Higgs donblets 9 might be indnced by radiative and finite-temperatnre effects, althongh it 
vanishes at the tree level. Withont any complex parameter, the phase 6 could be induced 
by loop effects of SUSY particles. At zero temperature, this idea of spontaneous CP 
violation was stndied in the minimal snpersymmetric standard model (MSSM) and it is 
pointed ont that there inevitably appears a psendoscalar boson as light as several GeV[9]. 
The same mechanism at finite temperatnres was snggested in Ref. [10] and extensively 
stndied by some of us[11]. They found that 6 conld be large only in the transient region 
between the symmetric and broken phase regions, with a psendoscalar Higgs whose mass 
satisfies ttia < 85GeV. This mechanism was appealing in that snch a large 6 can prodnce 
snfficient BAU and it does not indnce a large 6 in onr vacnnm, which is consistent with 
the bonnd from nentron electric dipole moment (nEDM). This scenario, however, is now 
excluded since the mass of the psendoscalar must satisfy uia > 90.1GeV[2]. 

Nonzero 6 is also indnced from complex parameters in the MSSM. Bnt magnitnde of 
explicit CP violation is constrained by nEDM measurements. For example, the physical 
CP phase relevant to the EDM must be as small as 0(10“^) when masses of the SUSY 
particles are weak scale, or they are heavier than ITeV when the phases are 0(1) [8, 12]. 
Recently it was observed that the spectrnm and interactions of the nentral Higgs bosons 
are affected by a large explicit CP violation in the third generation of sqnarks sector, 
which is not restricted by the nEDM constraints[13]. The anthors fonnd that when the 
imaginary part of the prodnct of fi and the Higgs-sqnark trilinear conpling is large, the 
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lightest Higgs boson iJi, which is composed of the scalar and pseudoscalar Higgs helds, 
becomes much lighter than the present bound llSGeV, but is hard to be observed, since 
its couplings to the gauge boson and to the bottom quarks is very small. 

One may think that this CP violation can generate sufficient BAU in the MSSM 
for parameter sets allowed by experiments. Our purpose is to investigate the effects of 
the explicit CP violation on the EWPT in the MSSM and to evaluate the magnitude 
of CP violation relevant to electroweak baryogenesis, which should be measured at the 
transition temperature. The organization of the paper is as follows. In Section 2, we 
formulate the effective potential of the Higgs helds including the one-loop corrections from 
the gauge bosons and the third generation of quarks and squarks, both at zero and hnite 
temperatures. The masses of the three neutral Higgs bosons and the charged Higgs boson 
are dehned as the derivatives of the effective potential at the zero-temperature vacuum. 
The mass formulas are almost the same as those in Ref. [13] except for inclusion of the 
gauge boson contributions. In Section 3, we study the EWPT for parameter sets, which 
are consistent with the mass bounds on the neutral lightest and charged Higgs bosons, 
in the absence of the explicit CP violation. Next we introduce the phase of the trilinear 
coupling, examine how the strength of the phase transition changes, and evaluate the 
magnitude of the CP violation relevant to electroweak baryogenesis. Section 4 is devoted 
to concluding remarks. We summarize the formulas for the Higgs masses in Appendices. 

2 Effective potential of the MSSM 

We consider the MSSM that has the following superpotential, 

W = £« {f'ABHiQADB - ffsKUB - . (2) 

Besides supersymmetric lagrangian, the low-energy MSSM contains the soft-SUSY-breaking 
terms 

^soft = - + h.c.) 

- - ml^B^ARdBR “ ml^BU^RUBR 

- e,, ^dQ^LdBR - ^ud^LUBR + h-c.] . (3) 

We calculate the effective potential of the Higgs helds by taking into account the one- 
loop contributions from the gauge bosons and the third generation of quarks and squarks. 
We consider the gauginos to be heavy enough to decouple so that the most dangerous 
contribution to nEDM from the gluino is negligible. The correction from the leptons, the 
other quark and squarks can be neglected because of their small Yukawa couplings to the 
Higgs helds. Now the ehective potential at zero temperature is given by 

T = 0) = Voi^d, $«) + AoVi^d, $.), (4) 

where $«) is the tree-level potential 

2 

( 5 ) 


Ho = + h-c.) 


gi + gf 


(4*4., - 4*4.)" + f 14*4. 
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and AoV^($d, $«) is the one-loop correction written as 


Aol" 


N, 


c 


327r2 


E 

q=t,b 


El 

i=i,2 


ml)' 


6471-2 


ml)' 


log 


ml 

M2 





( 6 ) 


Here m? , m? m| and are field-dependent masses defined by (19), (20), (21), 
(23) and (24), respectively M denotes the renormalization scale, which we choose such 
that the loop corrections vanish at the vacuum. The expression (6) is the same as the one- 
loop correction in [13], except for our inclusion of the gauge-boson contributions, which 
strengthen the first-order EWPT. 

It is well known that the masses of the Higgs bosons receive large corrections from the 
loops of the top quark and squarks[14]. Here, the masses of the Higgs bosons are defined 
by the second derivative of the effective potential at the vacuum. To evaluate them, we 
parameterize the Higgs fields by the vacuum (n^, Vu, 0) and fluctuation around it as 


$. = P2 


{vd + hd + ittd) 


^u = 


je 


1 

P2 




(7) 


In the following, we represent the quantities evaluated at the vacuum by (■■), that is, 
evaluated with all the fluctuations being set to zero. Requiring that the first derivatives 
of the effective potential with respect to the neutral Higgs fields evaluated at the vacuum 
vanish, that is. 


dhd / \ ohd / \ dad / \ dad / 


( 8 ) 


we can express and ml in (5), in terms of Re(m|e*®), tan/? = Vu/vd and the particle 
masses from the first equation, and in terms of the particle masses from the 

second equation. Now the mass-squared matrix of the neutral Higgs scalars is expressed 
as 


/ Ml MIp \ 

\ (Mlpf Ml J ’ 


(9) 


where 
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where each component is given in Appendix B. One can hnd that the pseudoscalar ele¬ 
ments are factorized as 


cot,). (-) 

so that, the unphysical Goldstone mode can be extracted by 

( (^d \ ^ ( cos(3 sin(3 W \ 

\ ttu J y — sin P cos P j y a j ' 


Hence, the mass-squared eigenvalues of the neutral Higgs bosons are the eigenvalues of 
the matrix 


M 


2 _ 
0 — 


(-M|)i2 

V ^ (•^|p)l2 


(-M|)i2 

(-^ 5)22 

(Af|p)2i 


^(Af|p)i2 ^ 

(-M|p)2i 
_1_ ( A/<2 A 


(13) 


In the presence of the CP violation which induces the scalar-pseudoscalar mixing, 
the pseudoscalar is no longer a mass eigenstate. In what follows, we use the mass of the 
charged Higgs boson as an input parameter, instead of the pseudoscalar. The mass matrix 
of the charged Higgs scalar has the form of 


/ (9^HeflF \ / tan/3 1 \ 

\dp-^dp-/\ 1 cotpj' 


Similarly, we can extract the Goldstone mode so that he mass of the charged scalar is 
given by 


mp± 


1 

sin P cos P 



(14) 


The detailed form of the mass is given in Appendix G, by which one can express Re(m|e*®) 
in terms of mp±. 

The effective potential at hnite temperature contains temperature-dependent correc¬ 
tions; 

Heff($d, T) = Kff($d, T = 0) + ApH(4>,, <h,; T), (15) 


where 


ArH = 


27r^ 


6/A!!t(hU3/y’"-w 


2^2 


2^2 
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Here the functions Ib{cP) and If{cP) are dehned by 

/b,f(o^) = r dx X? log (1 T . 


( 16 ) 


( 17 ) 
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The function /(a^) yields a^-term with negative coefficient when expanded for 1[15]. 
This qualitatively explains why the EWPT becomes first order with bosons whose held- 
dependent mass-squared behave as m^{v) ~ for small Because of its large Yukawa 
coupling, a stop with a vanishing soft mass makes the first-order EWPT stronger. In the 
following, we will not use the high-temperature expansion (m^/T^ ^ 1 ), but numerically 
calculate the integrations to study the EWPT quantitatively. 

An important constraint on the hnite-temperature effective potential comes from the 
requirement that the sphaleron process decouples immediately after the EWPT, in order 
for baryon asymmetry produced at the EWPT not to be washed out. If we denote the 
minimum of the effective potential as {vd,Vu,0), the first-order EWPT is characterized 
by the degenerate minima at the transition temperature Tq] {vc cos Pc,vc sin Pc,(^c) in 
the broken phase region and (0, 0, 6q) in the symmetric phase region. The sphaleron 
decoupling condition is now written as[16] 


Vc 


T, 


> 1 . 


c 


( 18 ) 


The difference between Oq and 6q is crucial to determine total amount of baryon asym¬ 
metry produced at the EWPT. The prohle of the bubble wall created at the first-order 
EWPT is derived from the equation of motion for the gauge-Higgs system with the effec¬ 
tive potential at Tc[17]. Then the boundary conditions for the Higgs helds are provided by 
{vc cos Pci Vc sin Pc, Oc) and (0, 0, 6 ^ 0 ). Since the bubble wall prohle smoothly interpolates 
between the degenerate minima, the phase 9 varies from Oq to dc at the phase boundary. 
If there is no local minimum of the effective potential at Tc which leads to the transi¬ 
tional CP violation[ 11 ], we expect that the phase monotonously varies so that the baryon 
number generated is characterized by Oc — 9 q. Once we choose the phase convention such 
that 6 = 0 at the zero-temperature vacuum, both Oc and Oq are dehnitely calculated from 
the effective potential in the presence of the explicit CP violation in the squark sector. 
As seen from ( 8 ), Im(m|) becomes nonzero and Oq = —6 with S = Arg(m|)^. Then Oc + S 
is another important quantity we must evaluate. As shown in various works, when Oc + S 
is 0{1), sufficient BAU is generated by the charge transport mechanism[18]. 


3 Numerical results 

There are many parameters in the model, some of which can be fixed by requiring the 
vacuum at zero temperature to be the prescribed one characterized by Vq = 246GeV, 
tan/? and 0. In practice, we determine ml and by use of the tadpole condition ( 8 ). 
Re(m 3 e*^) can also be determined, once the charged Higgs mass mH± is given by (44). 
In order to evaluate the right-hand-sides of these equations, one must prepare p and the 
parameters in the squark sector; soft masses and trilinear couplings At and 

Ah. Since the CP symmetry is violated by nonzero Im(pAi) and Im(pAfe), we take p 
to be real and regard the phases of the A-parameters as inputs. The EWPT is strongly 

"‘in the second paper of [17], the relation 6q = —6 is proved by use of the kink ansatz for the wall 
prohle. One can show that this relation generally holds, by use of the asymptotic expansion of the wall 
prohle in the symmetric phase. 
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first order when the soft mass of squark is very small, as noted in the last section. One 
usually choose ~ 0 and niq = C>(100)GeV to avoid too large corrections to the p- 
parameter. For dehniteness, we take = lOGeV, = lOOGeV, and several values 
of rriq larger than ITeV. As seen from the scalar-pseudoscalar mixing elements of the 
mass-squared matrix (13), the CP violation induced in the Higgs sector is enhanced for 
a larger Im(/iAJ. In the following, we take large values for p and |At|, but to avoid a 
color-and-charge-breaking vacuum, we keep pcotp = A = \At\ = Thus all the 

input parameters are tan/?, p, mH±, rriq and Sa = ArgA* = ArgA^. 

First of all, we turn off GP violation and calculate masses of the neutral Higgs bosons 
for various tan/?, mH±, p and rriq. Among parameter sets consistent with the present 
bounds on the lightest Higgs boson mass ttihi, we pick up several ones, for which the 
EWPT is investigated. We numerically calculate the effective potential (15), and search 
for the minimum by use of the downhill simplex algorithm. We dehne Tc as the tem¬ 
perature at which this minimum degenerates with the symmetric phase, and evaluate vc, 
tan Pc and Oq. Next we gradually increase Sa and examine Tq, vq and Oq + S. 

Before showing the numerical results, we roughly describe how the spectrum, CP 
violation and the strength of the EWPT depend on the parameters. If we increase mH±, 
which implies larger Re(m|), the masses of the neutral Higgs boson grow and scalar- 
pseudoscalar mixing decreases, since the diagonal elements of (13) increases. Then, as is 
well known, the EWPT becomes stronger for larger mH± ■ If we decrease tan P, which 
implies larger top Yukawa coupling yt and larger A for a hxed p, the EWPT becomes 
stronger and Higgs-mixing is enhanced. As for the effects of 6a, we expect that the 
strength of the first-order EWPT will be weakened. This is because nonzero 6a modifies 
n-dependence of the mass-squared of the lighter stop, which is roughly proportional to 
when = 0 and 6a = 0, as seen from (26). 

Now we present numerical results. In the absence of CP violation, the bounds on 
the mass of the lightest neutral Higgs excludes the most portion of theoretically allowed 
region[2]. Por small tan/?, the lower bound is the same as the MSM, rrih > 115GeV. 
Eor tan/? between 8 and 40, the lightest Higgs boson can be as light as 92GeV. Very 
large tan /^-region is excluded by the fact that GDF at Fermilab have not observed 66-pair 
from the lightest Higgs boson. As for the charged scalar, the bounds are satished with 
mH± > 90GeV for tan/? = 1 — 50. First of all, we turned off the CP violation and 
calculated the masses of the Higgs bosons for tan/? = 5, 10, 20 and 30. For each tan/?, we 
studied the EWPT at several points (p, mp) which are allowed by the Higgs mass bounds. 
An example of contour plots of the lightest Higgs mass in (m//±,m^)-plane is shown in 
Fig. 1. There is no parameter set with tan/? = 5 for which the EWPT is strongly first 
order satisfying (18). For tan/? = 10, 20, we found several parameter sets with small rriq 
for which the EWPT is strongly hrst order. The lightest Higgs boson mass is smaller than 
about llOGeV for such parameter sets. 

Next, introducing CP violating phase 6a, we studied the strength of the EWPT and 
the CP violation relevant to electroweak baryogenesis at several parameter sets for which 
the Higgs mass bounds are satished and the EWPT is strongly hrst order in the absence 
of the CP violation. In our convention, 6a — 7^ is disfavored by the b ^ s'y constraint [19]. 

®In practice, we also studied the case with /icot/3 A A, but the results are not so altered as long as 
the squarks do not have nonzero expectation value. 
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Figure 1: The lightest Higgs mass for tan/? = 10 and /i = ISOOGeV. The dot stands 
for parameter set for which the EWPT is strongly Erst order, and the cross for that not 
satisfying (18). 


For example, the results are shown in Fig. 2 for the point in Fig. 1 with mH± = ISOGeV 
and rriq = ITeV. Although we do not show the mass of the second-lightest Higgs boson, 

GeV rad 




h^ldeg) h^ldeg) 

Figure 2: In the left-hand figure, the (5A-dependences of vc (solid curve), Tc (dashed curve) 
and the lightest Higgs mass (dotted curve) are plotted for tan /? = 10, p = ISOOGeV, 
mH± = ISOGeV and mq = ITeV. In the left-hand figure, the dashed curve stands for 
5 = Aigiml), the dotted curve for 9c and the solid curve for 9c + 5 for the same parameter 
set. 

the mixing in the Higgs bosons becomes maximal for 5a > OOdeg. The strength of the 
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EWPT becomes too weak to satisfy the sphaleron decoupling condition (18) for Sa > 
40deg. The magnitude of CP phase relevant to baryogenesis is sufficiently large for 
Sa > lOdeg, in spite of small Higgs mixing. For a larger charged Higgs mass, the EWPT 
becomes stronger, while the Higgs mixing becomes smaller, as expected. The results 
for mH± = 200GeV are shown in Fig. 3. The strongly first-order EWPT persists for 
Sa < 55deg, while the magnitude of Oq + S decreases. We also examined the EWPT and 


GeV rad 



0 1 0 20 30 40 50 60 70 80 90 0 1 0 20 30 40 50 60 70 80 90 

5A(deg) 5A(deg) 

Figure 3: The same as Fig. 2 but with mH± = 200GeV. 

the CP phase for larger and found that the strongly Erst order EWPT persists for 

larger Sa, but \dc + 5| becomes very small. For example, the maximal value of \9c + 5| is 
about 0.02 for mH± = ITeV and 0.005 for mH± = 2TeV. 

For larger /i, the effect of Sa is expected to become stronger, that is, the EWPT is 
weakened for smaller Sa- The results for jj, = 2500GeV and rriq = llOOGeV are depicted 
in Fig. 4. The strongly first-order EWPT persists for Sa < 30degree, while \9c + S\ is 
0{1) for Sa — 20degree. If we increase mH±, the EWPT gets stronger so that (18) is 



6A(deg) 6A(deg) 


Fignre 4: The same as Fig. 3 bnt with /i = 2500GeV and rriq = llOOGeV. 
satisfied np to larger Sa, bnt the magnitnde of CP violation \9c + S\ decreases, as the 
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case with smaller /x. As an example of a larger tan/?, we show the results for tan/? = 20, 
fi = 2500GeV and rriq = 1220GeV in Fig. 5. As noted above, a larger tan/? implies 
a smaller top Yukawa coupling, so that the effect of CP violation in the stop sector 
decreases. In fact. Fig. 5 shows that the strongly first-order EWPT persists for a larger 
5a and \6c + 5\ is smaller than than the case with tan/? = 10 (Fig. 4). We also explored 


GeV 

125 n - 1 
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115 - 
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105 - 
100 - 
95 - 

90 -|-^- 1 -^- 1 -^- 1 -^- 1 - 

0 1 0 20 30 40 50 60 70 80 90 

5/i(deg) 


rad 



5A(deg) 


Figure 5: The same as Fig. 3 but with tan/? = 20, /i = 2500GeV and rriq = 1220GeV. 


other parameter sets with the same tan (3 and observed the same tendency as the case 
with tan /? = 10. 


4 Conclusion 

We have studied the EWPT and CP violation in the Higgs sector of the MSSM. For 
parameter sets consistent with the present bounds on the masses of neutral and charged 
Higgs bosons, we found strongly first-order EWPT when the lighter stop mass is less than 
that of the top quark and the lightest Higgs mass is less than about llOGeV for 8 < 
tan/? < 30, in the absence of CP violation. These results without CP violation are not 
new. Since the transitional CP violation cannot occur for parameter sets consistent with 
the updated Higgs boson mass bounds, viable CP violations for electroweak baryogenesis 
are among those in the complex parameters in the model. The relative phase between fj, 
and gaugino soft mass is essential for the scenario in which the charginos and neutralinos 
play the role of charge carriers, while the phase is strongly constrained by the neutron 
EDM bound. If the masses of charginos and neutralinos are found to be as heavy as ITeV, 
they cannot participate the baryogenesis, since they can hardly be excited at the EWPT 
temperature. Another source of CP violation is the relative phase between fi and the 
squark-Higgs trilinear coupling A. In particular, that phase in the third generation is free 
from the nEDM constraint, and is expected to play an important role in baryogenesis. 
Further, it was pointed out that the explicit CP violation in the stop sector can make 
the scalar-pseudoscalar mixing in the Higgs bosons very large so that the Higgs boson 
lighter than the present bound is allowed because of their couplings to the gauge bosons 
and 6-quarks[13]. We investigated how this phase affects the EWPT. We found that for a 
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larger phase, the scalar-pseudoscalar mixing increases, but the first-order EWPT becomes 
weaker. The EWPT cannot continue to be strong enough for successful electroweak 
baryogenesis before the mixing reaches the maximal. Hence the Higgs scalar lighter than 
the present bound, which is allowed with a large explicit CP violation, is not consistent 
with the strongly first-order EWPT which is essential for electroweak baryogenesis. We, 
however, found that the phase relevant for the baryogenesis can be sufficiently large for the 
parameter sets consistent with the present mass bounds of the Higgs bosons. Since baryon 
asymmetry produced by the electroweak baryogenesis depends on the phase difference of 
the Higgs sector between the broken and symmetric phase regions, which are separated 
by the bubble wall created at the hrst-order EWPT. In our phase convention, the phase is 
6c + 5. The phase can be 0(1) for < llOGeV, mH± < 200GeV and 8 < tan{3 < 20, 
and decreases for a larger mH±- 

The calculations here are done at the one-loop level. At hnite temperatures, the naive 
loop expansion is not always reliable, and the resummed perturbation with temperature- 
corrected masses will improve the approximation. In particular, the infrared behavior of 
the Higgs scalar loop must be treated carefully. This is the reason why the results based on 
perturbation in the MSM was modihed by the improved perturbation or nonperturbative 
lattice calculation. As for the MSSM, the EWPT is mainly controlled by the stop loops, 
and the one-loop results in Ref. [20] are consistent with the results obtained by improved 
perturbation[21] and nonperturbative lattice calculation[22, 23], for the parameters with 
small fi and A. We expect that the results in this work will not be altered by such 
improvements. 

The parameter region of the MSSM allowed by the Higgs mass bounds is now much 
narrower than that allowed theoretically. Within that region, strongly first-order EWPT 
is possible only for < nit and niHi < llOGeV, which corresponds to 8 < tan/? < 30. 
The explicit CP violation in the stop sector can induce large CP phase in the Higgs 
sector which is relevant to electroweak baryogenesis. Whether this scenario works or not 
depends on the spectrum of the lightest Higgs boson, charged Higgs boson and the lighter 
stop. The masses of these particles will be clarihed in the near future by LHG. If the 
lightest Higgs boson is heavier than llOGeV, the EWPT will be too weak to make the 
sphaleron process out of equilibrium. Even if the EWPT is strongly Erst order, the model 
with heavy charged Higgs boson and large fi requires another source of CP violation than 
that in the Higgs sector for successful electroweak baryogenesis. It might be the relative 
phases of p and the gaugino masses, which are constrained by nEDM experiments. For 
these phase to work, the masses of the charginos and neutralinos should be as light as the 
weak scale, and it will also be checked in the near future. 
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A Field-dependent masses 

In this section, we snmmarize the held-dependent masses of the quarks and squarks of 
the third generation, and that of the gauge bosons, which appear in the dehnition of the 
effective potential (4). These are functions of the neutral components of the Higgs fields, 
while the charged components are needed to calculate the mass of the charged Higgs 
boson. The masses of the quark of the third generation are given by 



Similarly, the top and bottom squark masses are 



The gauge bosons masses are 


hiz = {^(1 +hdf + huf + (23) 

= -^92 iVd + hdY + + hu)'^ . (24) 

The masses evaluated at the zero-temperature vacuum are given as follows: 
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where Vq = nj + v^, tan/A = Vu/vd- 

B Mass matrix of the neutral Higgs boson 


(27) 

(28) 


The calculation of the elements of the mass matrix is presented in Ref. [13]. We included 
the gauge-boson contributions to their calculation, so that we record their explicit forms. 
For later convenience, we introduce the following quantities: 
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Mf = ml — mL + cos(2/5), 
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Ffe = |/i|^ — Ff, cot/5, Qb = \Ab\^ — Rbta.\i(3. 

The elements of the mass matrix in the scalar sector are given by 
= Re(m 3 e*®) tan /5 -f m| cos^ /5 
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= Re(m 3 e*^) cot f3 + m| sin^ f3 


, ^c- I ^ 2 ■ 2 n 

+Ifa5 2 ’""=” 


-a;,M2 + 2||/*|"g. + 2 ||/,|'P. ^ 


log^ + 


+2m: 


2 /-a;iM/ + 2||/i| m- 2 , 

• --‘°s gg - 2 i*'‘i ‘“s m5 


|2 ^72 + 771 


^Wu + 2 ||/f|^PtCot/3 
2Am? 


, H , 2(91 + 9 ! 

' „ w?r m? 

log — k — mj log — 2 . _ 

tl ^ ^2 *2 6 j^^2 t 


2Aml 

0 


, xlvl + 2 \yhf Rkcotp 2, 2, A 2 

+-^T -2 - log TT^ - log TTT - ^mr 


-XtMf + 2 \ytf Qtf 


2(Am?)2 


-x,M^ + 2\y,\^P^f 
2(Am?)^ 


m! + m? m? 

-Ji -Aw^ 

2Am? m? 

r 12' 

m? + ml ml 

_A-A log—A 

2Ami ^ ml 
0 02 


I (,i + log ^ + 2 , 2 ^ log ^ 


log^ 

m~ 

b2 


2 2 

log —- + log —^ 
^ M2 ^ M2 


= —Re(m 3 e*^) — m| sin/9 cos/9 


a;fcM2 + ||/,|^(g,-P,) ^ 

Ami m? 

0 62 


+m^ cot f3 


XtM^ + 2\ytfR 


log —Y /9- 


-a;fcM2 + 2||/,|^Pfc^ ml 


+ ^m|sin/3cos/3 (\ytf- 


xlvdVu + 2\ytf Rt 


2Aml 


2Aml 

0 


2 9l + 9!\ i, H 


i2 772 + 771 


2 2 
r> m- ^ m~ 


ilLr 1 IL 7 

log —- + log —- 
^ M2 ^ M2 


xlvdVu + 2\ybf Rb 


log log 

' 2 , 2 , . 2 

log -j^ - log j;p - ^^b 


+VdV, 


{xtM^ + 2 \ybfR){-xtM^ + 2 \ytf Qt) 
2(Am?)2 


ml + ml ml 

1- tl -A 1q _h_ 

2 Ami ^ ml 

t 12 
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(xfcMfc + 2 ||/fe|^ Qb){-XbM^ + 2 \ybf Pb) 

2 {Amiy 


+VdVu 

z[ixmir 


m? + ml 



( 36 ) 


The matrix elements in the pseudoscalar sector are 
= (A^p)_tan/3, 


{■^p)22 ~ (-^ 
{Ml) 


12 

Pl^^cotp, 


= Re{mle^^) + 

12 IGtt^ [ Am? 


^ 9 

log T7i - log 


4mg cot /3 2 
Am? * 


1 - 


M 2 
ml + m 


2 


2 Am? 


, m, 
log 


t2 ' 


m; 


bfe 


Am? 

0 


7Tl~ 

Rb I ml log ^ - ml log 


M 2 


m^ 


62 


M2 


— Am? 


4mgtan/3 2 

Am? ^ 


m| + m| 

fei _^ 

2 Am? 

0 


log- 


m 


f>i 


m 


b2 , 


The scalar-pseudoscalar mixing elements are given by 
{mIp)^^ = tan/3, 

(A 1 |p )^2 = (-^|p)ni cot/3, 

(ai 


2 

SP 


21 
2 ; 


Nc j mfit cot2 /3 
87r2 


Am? 


^ 2+^1 


log 


m: 




m^ 

Am- 

0 


Ivb 


Mrfbog™ 


Am? 


m? -|- m? 

tl _^ 

2 Am? 


(37) 

(38) 


(39) 


log- 


mt 


m 


t2 ' 


m 


+ 2\yt\ Qb 


62 


Ami 

b 


ml + ml ml 

bi iQg . 61 


2 Am? 

b 


m 


^2 , 


{m 


SP 


^ Nc [ mllt 
21 87 r 2 \ Ami 


I |2 92 + 91 

ml - 3 — 


log- 


m: 


I -XtM^ + 2\yt\^ Qt 


m_ 


+ 


*2 


Amj 


^ ml + ml 


( 40 ) 

( 41 ) 


( 42 ) 


m~ 

_ u_ 

mj 

t2 ^ 


mlh tan^ /3 
Ami 

b 


9 I + 9 I 


log 


, -a;feM2 + 2||/b|V, 


m 


+ 


62 


Ami 

b 


ml + ml ml 


( 43 ) 


•which are all composed of the terms proportional to Im(/iAt)/Am? or Im(/iAfe)/Am? 


C Charged Higgs mass 

The calculation of the charged Higgs mass is a tedious task, and the method is described 
in Appendix of Ref.[13]. We present the result, which contain the contribution from the 
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gauge bosons: 

2 1 
mH± = 

+ 


Re(m3e*®) + 


sin (3 cos (3 

Nc \ 1 

sin P cos P \ Am? 


/Kv 


(m| - - m?J 

— {ii 12 ) 


+ \yt\^Rt m? log 


mj 


M2 


- 1 


Am? 




(m? -ml){ml - ml) 


+ \ybfRb ml log 


m 


fei 


M2 


{bi 


2 ||/t|/b| mtmb 
2 2 






m. 




- 1 


(44) 


where 

fim. 


- ) 
qk^ 

= ^ \ytybf VuVd 


2ml - Tv(M'i + Ml)mi + (mQ 


22 J 


( \yt? + \yb? - y ) {ml - {Mfjjim^ - (Mf) 


22' 


VuVd 


-{\yi? + \yi2?)\\yi? + \yA"-^-w\^^ 


2 92 


+ \yt\ \yt\ ^ 



Qk 


l)^ + M^{\y3-§ 



22 


+ \ytybf + 2m? + 2ml + 2ml - m^ - ||uo cos(2/3) 

+ ^ \ytAt\^VuVd i\yb? - y) {m\ - {MfjJ + ]- \ybAb\^VuVd i\yt? - y) {m\ - {MfjJ 


-- \ytybAtAb\^ VuVd - \ytyb\^ VuVd{RtRh + hh) 


—2mlml 


R - f) «< + 


1 . 

y \ytyb\ VuVdRe{AtAl) 


.P-|)«. 


2m| - 2m? - m^ + cos(2/3) 


( 45 ) 
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